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ABSTRACT 

We report here the detection of a weak magnetic field of 50-100 G on the 09.7 su- 
pergiant £ Ori A, using spectropolarimetric observations obtained with NARVAL at 
the 2m Telescope Bernard Lyot atop Pic du Midi (France). £ Ori A is the third O 
star known to host a magnetic field (along with 8 1 Ori C and HD 191612), and the 
first detection on a 'normal' rapidly-rotating O star. The magnetic field of £ Ori A 
is the weakest magnetic field ever detected on a massive star. The measured field is 
lower than the thermal equipartition limit (about 100 G). By fitting NLTE model 
atmospheres to our spectra, we determined that £ Ori A is a 40 M Q star with a radius 
of 25 R©and an age of about 5-6 Myr, showing no surface nitrogen enhancement and 
losing mass at a rate of about 2xl0 -6 M Q yr . 

The magnetic topology of £ Ori A is apparently more complex than a dipole and 
involves two main magnetic polarities located on both sides of the same hemisphere; 
our data also suggest that £ Ori A rotates in about 7.0 d and is about 40° away from 
pole-on to an Earth-based observer. Despite its weakness, the detected magnetic field 
significantly affects the wind structure; the corresponding Alfven radius is however 
very close to the surface, thus generating a different rotational modulation in wind 
lines than that reported on the two other known magnetic O stars. 

The rapid rotation of £ Ori A with respect to 9 1 Ori C appears as a surprise, both 
stars having similar unsigned magnetic fluxes (once rescaled to the same radius); it 
may suggest that the sub-equipartition field detected on £ Ori A is not a fossil remnant 
(as opposed to that of 9 1 Ori C and HD 191612), but the result of an exotic dynamo 
action produced through MHD instabilities. 

Key words: stars: magnetic fields - stars: winds - stars: rotation - stars: early type 
- stars: individual: £ Ori A - techniques: spectropolarimetry 



1 INTRODUCTION 

Stellar magnetic fields have been detected across a large 
range of spectral types. In solar-type and essentially all cool, 
low-mass (i.e., mid F and later) stars, magnetic fields (and 
activity) are observed, often featuring a complex topology, 
and are thought to be due to dynamo processes occuring 
within the outer convective layers. In hotter, more massive 

* Based on observations obtained at the Telescope Bernard Lyot 
(TBL), operated by the Institut National des Science de PUnivers 
of the Centre National de la Recherche Scientifique of France. 



stars with outer radiative zones, magnetic fields are also de- 
tected (with a significantly simpler topology though) but 
only in a small fraction of stars (e.g., the magnetic chem- 
ically peculiar stars among the A and late B stars). The 
situation might be similar (though less well studied) among 
early B and O st ars, with only two O stars (namely 6 1 Ori C 
and HD 191612. [Donati et alJlioO^ . l200rj ) and less than a 
handfull of early B-type stars (e.g., r Sco, (3 Cep, £ Cas, 
iDonati et aflbOOrj , |200ll : iNeiner et ail 120031 ) yet known as 
magnetic. 

Magnetic fields are nonetheless expected to play a 
significant role throughout the evolution of hot massive 
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stars, by modifying the internal rotation, enhancing chem- 
ical transport a nd mixing, and produc i ng enhance d sur- 
face abundances (jMaeder fc Mevnetl20oil2004l2005h . Mag- 
netic fields can al so dramatically influence t he way winds 
are launched (e.g.. lud-Doula fc" Qwockill2003) and the later 
phases of evolution (e.g., the collapse, iHeeer et al1l2005l ): 
a large number of observational phenomena (e.g., non- 
thermal radio emission, anomalous X-ray spectra, abun- 
dance anomalies, and Ha modulation) can also be explained 
(qualitatively at least) by the existence of a weak mag- 
netic field. Yet, the origin of magnetism in massive stars 
is still an open question, with a lively debate between 
two classes of models. While some models assert that dy- 
namo processes (either located in the convective core, e.g., 
ICharb onneau fc MacGregorll200ll. or acting w ithin the ra- 
diative zone, e.g., iMullan fc MacDon ald 2005) can produce 
the observed magn etic fields, some others claim that the fiel d 
is fossil in nature (|Ferrario fc Wickramasinghell2005l , [2006), 
being advected and amplified through the initial protostellar 
collapse. 

The limited knowledge that we have about the exis- 
tence and statistical properties of magnetic fields in massive 
O stars is mostly due to the fact that these fields are difficult 
to detect. Absorption lines of O stars are both relatively few 
in number in the optical domain, and generally rather broad 
(because of rotation or to so me other type of as yet unknown 
macroscopic mechanism, e.g.. lHowarth et al.lll997l ), decreas- 
ing dramatically the size of the Zeeman signatures that their 
putative fields can induce. The results obtained so far (on 
two stars only) suggest that magnetic O-type stars may be 
(i) slow rotators and (ii) may exhibit a peculiar spectrum 
with very regular temporal modulation. While this view may 
partly reflect an observational bias (magnetic detections be- 
ing easier on slow rotators) or a selection effect (observations 
often concentrating on peculiar stars first), null results re- 
cently reported on intermediate and fast rotators argue that 
this effect may be real. This question is nevertheless a key 
point for clarifying both the origin and evolutionary impact 
of magnetic fields in massive stars and therefore deserves 
being studied with great care. 

With the advent of the new generation spectropolarime- 
ters, such as ESPaDOnS at the Canada-France-Hawaii Tele- 
scope (CFHT) in Hawaii and NARVAL on the Telescope 
Bernard Lyot (TBL) in southern France, studies of stellar 
magnetic fields have undergone a big surge of activity; in 
particular, detecting magnetic fields of massive O stars (or 
providing upper limits of no more than a few tens of G) is 
now within reach. In this context, we recently initiated a 
search for magnetic fields in a limited number of 'normal' O 
stars, using NARVAL. 

One of our targets is £ Ori A, a 09.7 lb supergiant 
l|Maiz-Apellaniz et al.l[20o3 ) and the brightest O star at op- 
tical wavelengths. Evidence for azimuthal wind structura- 
tion (with an modulation timescale of about 6 d, com- 
patible with the rotation p eriod) is reported from both 
UV and optical lines (e.g., iKaper et al.l 1 19961 . Il999h and 
possibly due to the presence of a weak magnetic field. 
£ Ori A is also well-kno wn for its prominent X- ray emission, 
logLx/Lboi = -6.74 jBerghoefer et all 119971 ). The origin 
of this X-ray emiss ion is however still controversial; while 
ICohen et al.l (|2006h suggest that it is due to the classical 
wind-shock mechanism (with X-rays originating from cool- 



ing shocks in the acceleration zone), iRaassen et al. l (|2008f) 
invoke a collisional ionization equilibrium model and lPollockl 
l|2007l ) argue for collisionless shocks controlled by magnetic 
fields in the wind terminal velocity regime. For all these 
reasons, £ Ori A is an obvious candidate for our magnetic 
exploration program. 

In this paper, we report our spectropolarimetric obser- 
vations of £ Ori A and present the Zeeman detections we ob- 
tained (Sec.[2J). From the collected spectra, we reexamine the 
fundamental parameters of £ Ori A and discuss the observed 
rotational modulation to attempt pinning down the rotation 
period (Sec. [5J. We then carry out a complete modeling 
of the detected Zeeman signatures and describe the recon- 
structed magnetic topology (Sec. We finally summarise 
our results, discuss their implications for our understanding 
of massive magnetic stars, and suggest new observations to 
confirm and expand our conclusions (Sec. (5J. 



2 OBSERVATIONS 

Spectropolarimetric observations of £ Ori A were collected 
with NARVAL at TBL in 2007 October, as part of a 10-night 
run aimed at investigating the magnetic fields of hot stars. 
£ Ori A was observed during seven nights; altogether, 292 
circular-polarization sequences, each consisting of four indi- 
vidual subexposures taken in different polarimeter configu- 
rations, were obtained. From each set of four subexposures 
w e derive a mean Stok es V spectrum following the procedure 
of lDonati et al.l |l997j) , ensuring in particular that all spuri- 
ous signatures are removed at first order. Null polarisation 
spectra (labelled N) are calculated by combining the four 
subexposures in such a way that polarization cancels out, 
allowing u s to check that no spurious signals are present in 
the data IjDonati et al.lll997l . see for more details on how 
N is d efined). All frames were processed using Libre ES- 
pRIT (jDonati et all 1 19971 ; Donati et al., in prep.), a fully 
automatic reduction package installed at TBL for optimal 
extraction of NARVAL spectra. The peak signal-to-noise ra- 
tios per 2.6 kms -1 velocity bin range from 800 to 1500, 
depending mostly on weather condition s (see Tabl e [T]). 

Least-Squares Deconvolution (LSD; lDonati et al.lll997h 
was applied to all observations. The line list was constructed 
manually to include the few moderate to strong absorption 
lines that are not (or only weakly) affected by the wind. 
The strong Balmer lines, all showing clear emission from the 
wind and/or circumstellar environment at the time of our 
observations, were also excluded from the list. The C IV lines 
at 580.13 and 581.20 nm are used as reference photospheric 
lines from which we obtain the average radial velocity of 
£ Ori A (about 45 kms" 1 ); a few unblended absorption lines 
that are not blueshifted with respect to the reference frame 
by more than 15 kms -1 are also included in the list. We 
end up with a list of only six lines, whose characteristics are 
summarised in Table [2] 

From those lines we produced a mean circular polar- 
ization profile (LSD Stokes V profile), a mean check (N for 
null) profile and a mean unpolarized profile (LSD Stokes / 
profile) for each spectrum. All LSD profiles were produced 
on a spectral grid with a velocity bin of 7.2 kms -1 . Aver- 
aging together all LSD profiles recorded on each night of 
the seven nights of observation (with weights proportional 
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Table 1. Journal of observations. Columns 1—5 list the date, the range of heliocentric Julian dates, the range of UT times, the number 
of sequences and the exposure time per individual sequence subexposure, and the range of peak signal to noise ratio (per 2.6 kms -1 
velocity bin), for each night of observation. Column 6 lists the rms noise level (relative to the unpolarized continuum level and per 
7.2 kms" 1 velocity bin) in the circular polarization profile produced by Least-Squares Deconvolution once averaged over the whole night 
(Section [2}. The rotation cycle (using the ephemeris given by Eq. [TJ is listed in column 7. 



Date HJD UT t cxp S/N <r LSD Phase 

(2007) (2,454,000+) (h:m:s) (s) (lQ- 4 / c ) 



Oct. 


18 


391.53928- 


-391.69519 


00:52:56- 


-04:37:25 


48 x 4 x 20 


780- 


- 990 


0.41 


0.648- 


-0.670 


Oct. 


19 


392.69363- 


-392.72343 


04:35:04- 


-05:17:59 


8 x 4 x 40 


1010- 


-1080 


0.87 


0.813- 


-0.817 


Oct. 


20 


393.54410- 


-393.72674 


00:59:39- 


-05:22:38 


44 x 4 x 40 


1220- 


-1470 


0.28 


0.935- 


-0.961 


Oct. 


21 


394.46543- 


-394.66510 


23:06:15- 


-03:53:46 


48 x 4 x 40 


810- 


1460 


0.28 


1.066- 


-1.095 


Oct. 


22 


395.49302- 


-395.69959 


23:45:52- 


-04:43:19 


48 x 4 x 40 


1090- 


-1480 


0.27 


1.213- 


-1.242 


Oct. 


24 


397.47086- 


-397.67036 


23:13:45- 


-04:01:01 


48 x 4 x 40 


1030- 


-1480 


0.27 


1.495- 


-1.524 


Oct. 


25 


398.50205- 


-398.70310 


23:58:34- 


-04:48:03 


48 x 4 x 40 


1200- 


-1470 


0.27 


1.643- 


-1.671 



Table 2. Lines used for Least-Squares Deconvolution. The line 
depths (column 3) were directly measured from our spectra while 
the Lande factors (column 4) were derived assuming LS coupling. 



Wavelength Element 
(nm) 



Depth Lande 
(/ c ) factor 



492.1931 


He I 


0.44 


1.000 


501.5678 


He I 


0.37 


1.000 


541.1516 


He II 


0.35 


1.000 


559.2252 


O in 


0.40 


1.000 


580.1313 


C iv 


0.20 


1.167 


581.1970 


C iv 


0.15 


1.333 
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Figure 1. LSD Stokes V (top), null N (middle) and Stokes I 
profiles of £ Ori A on 2007 Oct 24. The V and N profiles are 
expanded by a factor of 500 and shifted upwards by 1.2 and 1.1 
for display purposes. A clear Zeeman signature is detected in the 
red line wing while the null profile shows no signal. 



to the inverse variance of each profile) yields relative noise 
levels of 0.27 (in units of 10 _4 7 c ) except on the first two 
nights (where the noise reaches 0.41 and 0.87). On Oct. 24, 
the detection probability exceeds 99%, with a reduced-x 2 
value (compared to a null-field, V — profile) of 1.33; the 
corresponding Stokes V (and null N) LSD profiles are shown 
in Fig. [T] Similar (though less clear) Zeeman signatures are 
also observed during the other nights. 



3.1 



PARAMETERS AND ROTATION OF < Ori A 
The photosphere and wind of £ Ori A 



We reexamine the spectral properties and fundamental pa- 
rameters of ( Ori A, using the new recorded spectra. 

We have performed the spectrum analysis using model 
atmos pheres calculated with the unified model code CMF- 
GEN (|Hillier fc Miller! Il998f l. CMFGEN provides a consis- 
tent treatment of the photosphere and the wind, thus of- 
fering useful insights into the wind properties while provid- 
ing a realistic treatment of photospheric metal-line blanket- 
ing. The code solves for the atmospheric structure, non-LTE 
populations and the radiation field, in the comoving frame 
of the fluid. The fundamental stellar parameters (T e g, logg, 
Ri, and abundances) must be specified at this step, together 
with the mass-loss rate and velocity law. After convergence 
of the model, a formal solution of the radiative transfer equa- 
tion is computed in the observer's frame, thus providing the 
synthetic spectrum for comparis on to observations. For more 
details on this code, we refer to iHillier fc Milieil (|l998l ) and 
iHillier et all (|2003l ). CMFGEN does not solve the full hy- 
drodynamics, but rather assumes a density structure. We 
use a hydrostatic densi t y structure computed with TLUSTY 
i|Hubenv fc Lanzll 19951 : fcanz fe Hubeny||2003h in the deeper 
layers, while the wind regime is described with a standard 
/3-velocity law. The photosphere and the wind are connected 
below the sonic point at a wind velocity of about 15 kms -1 . 

Radiatively driven winds are intrinsically subject to in- 
stabilities, resulting in the formation of discrete structures 
called "clumps". Both observational evidence and theoreti- 
cal arguments foster the concept of highly-structured w inds 
(|Eversberg et all Il998l ; iDessart fc Owockil 120031 , l2005h . To 
investigate spectral signatures of clumping in the wind of 
£ Ori A and its consequences on the derived wind parame- 
ters, we have constructed clumped wind models with CMF- 
GEN. A simple, parametric treatment of wind clumping is 
implemented in CMFGEN, which is expressed by a volume 
filling factor, /. It assumes a void interclump medium and 
the clumps to be small compared to the photons mean free 
path. Clumps start t o form in the wind at velocities higher 
than v c i. We refer to IHillier et all ((2003) for a detailed de- 
scription of wind-clumping. 

We adopted a value of the clumping filling factor of 
/ = 0.1. Test models using / = 1 revealed little change of 
the Ha profile. As for v c i , models with values from 30 (see 
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iBouret et al.ll2005h to 400 kms -1 showed a larger and larger 
shift of the central absorption component of Ha towards 
shorter wavelengths. The best match was obtained for v c \ ~ 
200 kms" 1 . 

A depth-independent microturbulent velocity is in- 
cluded in the computation of the atmospheric structure (i.e., 
temperature structure & population of individual levels). 
We chose a value of 5 kms -1 as the default value (see 
iMartins et al.| [2002). For the computation of the detailed 
spectrum resulting from a formal solution of the radiative 
transfer equation (i.e., with the populations kept fixed), a 
depth dependent microturbulent velocity was adopted. In 
that case, the microturbulent velocity follows the relation 
v t urb(r) = v m in + (?Wx — Wmin) v{r)/voc where w min and u max 
are the minimum and maximum microturbulent velocities, 
and Voa the terminal wind velocity. In this formulation, v m m 
is technically equivalent to £t, the microturbulence in the 
photosphere. We considered several values of « m in, searching 
for consistent fits for the photospheric lines. This is obtained 
for Vmi n — 10 kms" 1 in the photosphere. u max = 0.1 Uoo was 
adopted at the top of the atmosphere. 

The effective temperature T e s was derived from the ra- 
tio of He II to He I lines as usually done for O stars. Com- 
paring with models from 27,500 to 31,500 K with 1,000 K 
steps, we find that T e g = 29, 500 K provides the best fit (see 
Fig- El- Given the high quality of the observed spectrum, we 
can clearly exclude the 27,500 and 31,500 K models. Models 
at 28,500 and 30,500 K already show significant deviations 
with respect to the best fit. We thus (conservatively) con- 
clude that our estimate is accurate to ±1,000K. The lumi- 
nosity was derived from the observe d magnitude. C OH A 
has my = 1.76 and bib = 1.59 (|Mafz-Apellaniz et al.l 
12004!) . which for an in trinsinc color index (B-V)o of —0.26 
i Martins fc Plej|2006h . corresponds to an extinction Ay = 
0.28 mag (assuming Ry=3.1). We assumed that the dis- 
tance to Orion is equal to d = 414 ± 50 pc (the uncertainty 
being taken as the dispersion among th e various recent mea - 
surements published in the literature. iMenten et alj|2007h . 
The luminosity of £ Ori A can thus be estimated from my, 
Ay, d and the b olometric correction ( equal to —2.73 for 
T efi = 29,500 K. IMartins fc Pled [20061 ). We end up with 
L = 10 5 ' 64 Lq. Combining the uncertainty in T c ffand d leads 
to an uncertainty of about 0.15 dex for L. The corresponding 
radius is thus R — 25 ± 5 R© . We also derived log g from the 
shape of H/3, H<5 and He. We computed models with \ogg 
ranging from 3.0 to 3.75 with 0.25 dex steps. For the three 
lines, the best match is obtained for \ogg — 3.25. The quan- 
tification of the goodness of the fit by means of \ 2 indicates 
that the uncertainty is of about 0.1 dex. From the estimates 
of log g and R, one gets M = 40 ± 20 M Q . 

In Fig. |3J we show the position of £ Ori A in the HR 
diagram. The two other known magnetic O stars (8 1 Ori C 
and HD 191612) are also reported for comparison. We de- 
rive an age of 5- 6 Myr for C Ori A from the Geneva evo- 
lutionary tracks l|Mevnet fc Maederl 120031 ) . A simple lin- 
ear interpolation between the evolutionary tracks gives a 
mass of 39±8 Mq, in good agreement with (and more ac- 
curate than) the spectroscopic mass derived above. At first 
glance, £ Ori A thus appears as an evolved couterpart of 
both 6 1 Ori C and HD 191612. 

Abundances of a few elements can also be derived. Fig. [4] 
shows such determinations for CNO, i.e., the most important 
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Figure 3. HR diagram with the position of £ Ori A, 9 1 Ori C and 
HD 191612 indica t ed. P aram eters for the las t two s tars are from 
ISimon-Di'az et alj 1120061 ) and lWalborn et al.l i2003t) re spectively. 
Evolutionary tracks are from lMevnet fc Maederl l|2003h . 



elements to constrain stellar evolution models. These mod- 
els predict mixing of CNO-processed material to be more 
efficient in fast-rotating stars; accordingly, these stars reveal 
larger CNO s urface abundance anomali es at an earlier stage. 
For instance, iMevnet fc Maederl (|2003l ) predict for a solar- 
composition, 40 Mq star with an age of 5 Myr that CNO 
surface abundances should be equal respectively to 0.5, 5 
and 0.6 times their original values. Fig. |4]clearly shows that 
we do not observe the large N enrichment. The observed 
N III lines are best matched with an abundance close t o the 
value found in the Orion nebula llEsteban et al. 2004) and 



Jue found in the Orion nebula (Estcban ct al. 2004) anc 
Orion main-sequence B stars (|Cunha fc Lambertlll994l ) 
Fitting C in lines gives uneven results: while some lines are 
better fitted with the Orion nebula abundance, others indi- 
cate a small depletion. Finally, most (but not all) O II lines 
indicate that O is underabundant by a factor ~2. The solar 
composition assumed for all other elements give good fits 
to the observed spectrum (e.g., the 448.0 nm Mg II and the 
various Si ill and Si IV lines). 

We used archival IUE spectra to measure the wind 
terminal velocity from the blueward extension of the strong 
UV P-Cygni profiles and found v x = 2,100 kms" 1 . To 
estimate the mass loss rate, we relied on Ha only. As Ha 
is varying with time (see below), we tried to fit the two 
profiles respectively featuring the strongest and weakest 
emission in the far wings (i.e., the part of the profile that 
can be most reliably fitted with our ID wind model) to 
derive the range of M over the rotation cycle, yielding 
values of M = 1.4 - 1.9 x 10" 6 M yr" 1 (see Fig. 01. We 
achieved a better match to Ha adopting a fast velocity law 
with /3 = 0.8. Fitting the position of Ha central absorption 
requires that clumping starts rather high in the wind (at 
velocities of about 200 km/s). This is larger than what 
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Figure 2. Modeling the 447 nm He I and the 454 nm He II lines with atmospheric models corresponding to different T e fj values (2 left 
panels; logg = 3.25), and the H<5, H/3 and He lines with models corresponding to different log g values (3 right panels; T e g = 29,500 K). 
For this plot, we built a mean spectrum out of the spectra obtained on the night of October 18, 2007. 



iBouret et all <|2005l l found for early O supergiants. A 
summary of our results from this spectroscopic analysis is 
presented in Table [3] 



such that the overall Ha profile reaches the peak intensity of 
the strongest observed Ha emission (corresponding to phase 
1.082). We find that M < 2.5 x 10~ 6 Moyr" 1 (note that 



Previous determinations of the physical p ro pertie s 
of £ Ori A were made by lLamers fc Leithererl (| 1993T ) . 
They gave T cff 
actually based on 



30, 900 K, but this estimate was 
the effective temperature scale of 



IChlebowski fc Garmanvl |l99ll) and not from a direct anal- 
ysis of the star with atmospheric models^; our finding that 
T c ff = 29, 500 K can thus be regarded as a sign i ficant im- 
provement with respect to lLamers fc Leithererl (|l993h . In 
the same study, the luminosity of £ Ori A is estimated to 
L — 10 5 ' 9 1/0, as a result of the larger distance they assumed 
(500 pc) and to a larger T e g (and thus a larger bolometric 

correction). Again, our estimate is more robust. 

Concerning the abundance patterns, iRaassen et alj 

(2008) recently derived solar C and N content as well as a 
small O depletion from X-ray spectra. This is consistent with 
our measurements for these elements. In particular, their 
data reveal the absence of strong N enrichment. They also re- 
port minor Mg and Si enrichment, but this is not confirmed 
by our results. Finally. lLamers fc Leithererl (|l993l ) derived a 
mass loss rate of 2.5 x 10~ 6 M0yr _1 from radio measure- 
ments; scaled to a distance of 414 pc (they assumed 500 pc), 
this corresponds to M — 1.9 x 10 -6 Mq yr _1 , in good agree- 
ment with our determination. We set an upper limit to the 
actual mass-loss rate of £ Ori A by searching for the model 



1 It is now well established that the T e ff scale of O stars has been 
revised downward, ( e.g., |MjjUns_et_al. |2002| ; ICrowther et alj 
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X[nm] 
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I2002L iRepolust et al j|2004 iMartins et al.ll2005l) 



Figure 5. Estimating mass loss from Ha profiles featuring weak- 
est (left) and strongest (right) emission in the far wings, over our 
7d campaign. Matching the broad profile wings yields mass loss 
rates in the range 1.4— 1.9xl0 — 6 Mq yr — 1 . The two weak lines in 
the red wing are from C II. 
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Figure 4. Determining the CN (top) and O (bottom) abundances of £ Ori A. In the top panels, the synthe tic spectrum correspond ing 
to CN abundances expected for a 40 Mg star at an age of 5 Myr (according to the evolutionary tracks of iMevn et fc Macdcr 200^) is 
shown in red. In the bottom panels, synthetic spectra with O II lines are shown for models with O = 0.5, 0.75 and 1.0 Oq (evolutionary 
models predicting O = 0.65 O0 at 5 Myr). For this plot, we built a mean spectrum out of the spectra obtained on the night of October 



Table 3. Summary of stellar properties of ( Ori A, including pho- 
tospheric and wind parameters derived from the modeling with 
CMFGEN. Abundances are expressed relative to Hydrogen. See 
text for a discussion of the photospheric abundance patterns rel- 
ative to the initial/local content. 



Spectral Type 


09.7 lb 


Distance [pej 


414. ±50 


± IU bCLli 1U11 ll_J*_l \_1 


7.0 ±0.5 


11 Q1 Tl 1 \\i\Y\ C 1 
U oil! L 1V111 o 


110. ±10 


Vmac [kms -1 ] 


93. ±9 


Inclination angle i [°] 


40. 


T cS [K] 


29,500 ±1000 


logg [cgs] 


3.25 ±0.1 
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C/H 
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N/H 


6.0 ± 1.8 x 10 -5 


O/H 
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the synthetic profile then strongly overestimates the wings 
strength). We refer to Sec. [5] for a discu ssion of the origin of 
the H a emission peak. Note finally that lLamers fc Leithereri 
l| 19931 ) value assumes no wind-clumping (/ = 1), while we 
adopt / = 0.1 for the optical study. 

3.2 Temporal variability and rotation 

Through the Fourier transform of the 580 and 581 nm C IV 
and O in 559 nm line profiles (averaged over all lines and 
nights), we can obtain an accurate estimate of the rotational 
broadening parameter v sin(i) (by matching the position of 
the first zero in the Fourier profile, when this first zero is 
visible), and of the additional macroscopic turbulent velocity 
broa dening lines (often far) beyond their rotation profiles 
(e.g. lGrav||l98lD . We find wsin(i) = 110 ± 10 kms -1 (see 
Fig. [6]l ; for the macroturbulence profile (assumed gaussian) , 
we find v mac — 93 ± 9 kms -1 (corresponding to a full width 
at half maximum 110 ± 10 kms -1 ). 

Given the estimated radius (25 R©) and rotational 
broadening (110 kms -1 ) of £ Ori A, we conclude that its 
maximum rotation period P TO t / sin i is e qual to 11.5 d. Sev - 
eral papers in the refereed literature (e.g lKaper et aLlll997l ) 
mention a possible rotation period (or half-rotation period) 
of 6 d for £ Ori A from variations of Balmer lines. Look- 
ing at how Balmer lines evolve with time during our own 
observations (see Fig. [7] first two left panels) also suggests 
rotation periods of 6 to 8 d depending on which portion of 
the line profile we focus on. The central absorption of H/3 
and Ha (both maximum on Oct 18 and Oct 24, i.e. cycles 
0.66 and 1.51) are apparently varying with a timescale of 
about 6 d, in agre ement with the published estimate (e.g., 
Kaper et al.lll997n . The red emission of Ha and the red 
wing of H/3 (minimum on Oct 18 and Oct 25, i.e., cycles 
0.66 and 1.66) suggest a slightly longer period. The far blue 
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f Ori, Fourier transform of C IV & III profiles 
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Figure 6. Fourier transform of the average 580 and 581 nm C IV 
and 559 nm O III line profiles of f Ori A. The full line is obtained 
from the observed profiles, the dashed line corresponds to the 
model. 

wing of both Balmer lines shows evidence of excess absorp- 
tion on Oct 19 and Oct 22 (i.e., cycles 0.82 and 1.23) but 
not on Oct 25 (cycle 1.66), suggesting a (half) period of 
about 3.5 d; the far red wing of both lines also show distinct 
excess absorption on Oct 20 (cycle 0.95). As a matter of 
fact, we find that variablity is present in a large majority of 
lines. Most lines exhibit a clear time variable blueshift, max- 
imum on Oct 19 and Oct 24 (cycles 0.8 and 1.5), i.e. on a 
timescale of at least 4-5d, as showed for the case of He I 492 
nm Fig. [7] (third panel from the left). The average blueshift 
(with respect to C IV 580-581 nm) is —16 kms -1 , while the 
peak-to-peak maximal amplitude is —17 kms -1 , well below 
the radi al velocity variations induced by the companion of 
C Ori A l|Hummel et alj|2000r i. 

We also detect temporal variations in the 569.59 nm 
C ill double-peak emission line (see Fig. [7] right panel) where 
the relative intensity of both peaks vary from one night to 
the next. In our observations, the red peak features maxi- 
mum emission on Oct 20 and Oct 23 (cycles 0.95 and 1.51) 
while the blue peak shows maximum emission on Oct 22 (cy- 
cle 1.23), apparently in antiphase with the red peak of the 
same line and in phase with the excess absorption episode 
occuring in the far blue wing of both Balmer lines. The in- 
tensity ratio of both peaks is varying on a timescale of about 
3-4 d and may potentially be a good indicator of the rota- 
tion half-period. 

Given that we detect clear modulation on a period of 
about 3.5 d both in the far blue wings of both Balmer lines 
and in emission peak intensity ratio of the C ill line, we 
interpret the ~7 d period as the rotatio n period (rather than 
half the period as iKaper et al.|[l997l ). The 4-5d timescale 
seen in most photospheric lines is potentially also compatible 
with a 7d rotation period (e.g. if the two maximum blueshifts 
are unevenly spaced in rotation phase). 

We therefore assume in the following a 7 d timescale to 
phase our data, and use the following ephemeris to compute 
rotational phases: 

HJD = 2454380.0 + 7.0E. (1) 
We further discuss below the determination of the rotation 
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Figure 7. Temporal variations of the H/3 (left), Ha (second from left), He I 492 nm (third from left) and 569.59 nm C III (right) 
lines throughout our observing run. The time-averaged profile is plotted in red to emphasise variations. The rotational cycle of each 
observation (assuming a rotation period of 7 d) is written next to each profile. All profiles are shown in the star's rest frame (i.e., shifted 
by 45 kms -1 with respect to the heliocentric rest frame). 



period, using the detected Zeeman signatures to probe rota- 
tional modulation. Assuming a rotation period of about 7 d 
implies that the star is seen at an inclination (with respect to 
the rotation axis) of about 40°, and that the equatorial ve- 
locity of £ Ori A is about 170 kms -1 . While this is not quite 
extreme rotation by O star standards, this is already much 
higher than the two magnetic stars known to date (8 1 Ori C 
and HD 191612, both featuring equatorial rotation velocities 
lower than 30 kms -1 ). 



4 MODELLING THE MAGNETIC TOPOLOGY 
OF ( Ori A 

To model the Z eeman signatures o f £ Ori A, we use the 
imaging code of iDonati et al. (2006). The magnetic topol- 
ogy at the surface of the star is reconstructed as a spherical- 
harmonic expansion, whose coefficients are adjusted (with 
a maximum-entropy image reconstruction code) to ensure 
that the synthetic Zeeman signatures corresponding to the 
reconstructed magnetic topology match the observed ones 
at noise level. The magnetic image that we derive can thus 
be regarded as the simplest topology compatible with the 
data. 

This new imaging method has several advan tages 
with respect to t he old er one of H rown et all ill99ll ) and 
IDonati fc Brownl (Il997l ). The reconstructed field is directly 
expressed as the sum of a poloidal and a toroidal field. More- 
over, we have a direct and obvious way of constraining the 
degree of complexity of the reconstructed field topology by 



limiting the spherical harmonic expansion at a given maxi- 
mum I value, depending on the quality and temporal sam- 
pling of the available Zeeman data. 

For this modeling, we assume that «sin(z) = 110 kms - 
and i = 40°, as derived in Sec. 13.21 We also assume that the 
star and its magnetic topology are rotating as a solid body 
with a rotation period of 7 d; different values of the rota- 
tion period are also used to evaluate how much the result 
we obtain is sensitive to this parameter (see below). The 
line profile model (including macroturbulence broadening) 
used to describe the synthetic Zeeman signatures is the same 
as that introduced in Sec. [3] to describe the observed pho- 
tospheric lines, i.e., a simple gaussian at an average wave- 
length of 500 nm and with full width at half maximum equal 
to 110 kms -1 . Zeeman signatures are obtained by assuming 
the weak field approximation and an average Lande factor 
of 1.1. 

The complete set of Zeeman signatures and their cor- 
responding null profiles are shown in Fig. [8] along with the 
maximum entropy fit to the data, assuming either a simple 
dipole field or a more complex magnetic geometry (limited 
to £ = 3). The reduced xt associated to the global set of 
V and N profiles with respect to a non magnetic (B = 0) 
model are equal to 1.25 and 0.99 respectively (for a total 
number of 392 data points), indicating that the magnetic 
signal in the V profiles is unambiguously detected at a 10cr 
level while no significant signal is observed in the N profiles. 

Assuming that the star hosts a simple dipole field pro- 
vides a better fit to the data; however, the resulting xt 
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Figure 8. Observed (black) and modeled (red & blue) Stokes V 
signatures (left) and null TV profiles (right) of £ Ori A. The blue 
line corresponds to a simple dipole, while the red one corresponds 
to a more complex field having I = 3. A clear Zeeman signal is 
detected and consistently modeled in the Stokes V profiles. The 
rotational cycle of each observation (assuming a rotation period 
of 7 d) is written next to each profile. A la error bar is also 
plotted left to each profile. 



(equal to 1.14) is still significantly larger than 1, indicat- 
ing that the magnetic topology of £ Ori A is likely more 
complex. Using a spherical harmonics series expanded up to 
£ = 3 provides a unit fit to the data, i.e., is successfull 
at reproducing the data down to noise level; in particular, it 
provides a much better fit to the data at cycle 1.51 on the 
red side (negative dip) of the line profile (see Fig.[SJ). Unsur- 
prisingly, carrying out the same analysis on the N spectra 
only yields flat synthetic profiles, the corresponding \v f° r 
a non-magnetic model being already below 1. 

The dipolar and £ = 3 magnetic topologies derived from 
the data are both shown in Fig. [9] The reconstructed dipole 
has a strength of 61 ± 10 G, is roughly perpendicular to 
the rotation axis (inclination angle /3 = 83 ± 10°) with the 
positive pole facing the observer at phase 0.42 ± 0.03. The 
second (more complex) magnetic topology shows more con- 
centrated features (where the field reaches as much 100 G) 
and is mainly poloidal (the toroidal component containing 
less than 5% of the reconstructed magnetic energy). Given 
the limited resolution we have access to on the star (about 
six resolution elements across the equator given the fairly 
large width of the local profile), the moderate accuracy to 
which the Zeeman signatures are detected and the moderate 
phase coverage of the collected data, there is no real point at 



carrying out reconstructions with even more complex field 
topologies. 

We carried out the same analysis for different values of 
the rotation period (and corresponding inclination), looking 
for the period that produces the magnetic image with the 
smallest information content and a = 1 fit to the data. 
From this process, we derive that 7.5 d is a marginally better 
rotation period, with a la error bar equal to about 0.5 d. 
The magnetic topologies derived for this value of the rota- 
tion period are very similar to those of Fig. O Using this cri- 
terion, we also find that periods of 12 to 14 d are less likely 
to be the true rotation period of £ Ori A. Note that this 
constraint on the rotation period depends on the assumed 
magnetic geometry (here limited to an £ < 3 spherical har- 
monic expansion) given the limited span of our observations 
(7 d). 



5 DISCUSSION 

We report in this paper the detection and the first modeling 
attempt of the weak large-scale magnetic field of the 09.7 
supergiant £ Ori A. We detect a field that corresponds to 
local surface magnetic fluxes of only a few tens of G. The 
field is everywhere lower than 100 G, making it (by far) 
the weakest magnetic field ever reported in a hot massive 
star (|Donati et al.lll99d ; lAuriere et al.| [2007). In particular, 
this magnetic field is weaker than the thermal equipartition 
limit, equal to about 100 G for £ Ori A; this is the first sub- 
equipartition field unambiguously detected in a hot star. The 
magnetic chemically peculiar stars all show fields large r than 
their thermal equipartition limit (lAuriere et al . 2007). This 
detection also brings the number of known ma gnetic O stars 
to th ree, with £ Ori A thus .joining 8 1 Or i C (|Donati et al.l 
120021 ) and HD 191612 (|Donati et al.l 120061 1 . This is also the 
first magnetic detection in a 'normal' rapidly rotating O 
star. 

The detailed spectral modeling of £ Ori A provides 
T eS = 29, 500 ± 1,000 K and log 5 = 3.25 ±0.10 with normal 
abundances. It follows that £ Ori A is a 40 Mq star with 
a radius equal to about 25 Rq, seen from the Earth at an 
inclination angle of 40°. With an age of about 5-6 Myr, 
£ Ori A essentially appears as an evolved counterpart of 
both 8 1 Ori C and HD 191612. Given its evolutionary stage, 
£ Ori A is expected to show significant N enrichment at its 
surface (as well as moderate C and O depletion); the normal 
nitrogen abundance that we measure is thus surprising. It is 
tempting to suggest that magnetic fields may play a role in 
this process; this is however not what th e first evolutionary 
mode ls including magnetic field predict l|Maeder fc Mevnetl 
2005). More work (both on the observational and theoretical 
side) are required to investigate this issue further. From a 
fit to Ha, we estimate that the mass-loss rate is about 1.4- 
1.9xl0 -6 Moyr -1 . From the temporal variability of spec- 
tral lines and the modulation of Zeeman signatures, we find 
that the period of £ Ori A is about 7d. This is compatible 
with the usin(i) that we measure (from the Fourier shape of 
the photospheric C IV lines) and the radius that we derive 
(from the spectral synthesis), provided the star is view at 
intermediate inclinations (i = 40 °). 

Given that £ Ori A is typically 3 and 1.4 times larger 
m size than 9 1 Ori C and HD 191612 respectively, we find 
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Figure 9. Reconstructed magnetic topology of £ Ori A, assuming cither a dipolc magnetic field (top) or a more complex field having 
£ = 3 (bottom). Only the second (lower) topology provides a = 1 fit to the Stokes V data. In both cases, the three field components 
are displayed from left to right (flux values labelled in G). The star is shown in flattened polar projection down to latitudes of —30°, 
with the equator depicted as a bold circle and parallels as dashed circles. Radial ticks around each plot indicate phases of observations. 

that its overall unsigned magnetic flux (i.e., the integral of 
the absolute value of the magnetic field over stellar sur- 
face) is slightly larger (by a factor of about 1.5) than that 
of 9 1 Ori C but much smaller than that of HD 191612 (by 
about an order of magnitude) . 

The extremely long rotation period of HD 191612 
(about 538 d) suggests that the magnetic field is likely re- 
sponsible for having dissipated (through confi ned mass loss) 
most of the angular momentum of HD 191612 (|Donati et al.l 
2006). The slow rotation rate and extreme youth of 1 Ori C 
also suggests that primordial magnetic fields pervading the 
parent molecular cloud must have a strong impact onto the 
angular momentum dissipation throughout the cloud col- 
lapse, in qualitative agreement with what numer ical simu- 
lations predict (e.g., iHennebelle fc From ane 200jl). In this 
context, one would expect £ Ori A to rotate, if not as slowly 
as HD 191612 (whose intrinsic magnetic flux is much higher), 
at least more slowly than Ori C (whose intrinsic mag- 
netic flux is similar) given its later evolution stage; this is 
however not what we observe. No more than speculations 
can be proposed at this stage. One possibility is that the 
magnetic field of £ Ori A is not of fossil origin (as opposed 
to that of 1 Ori C and HD 191612) but rather dynamo 
generated, making the rotational evolution of £ Ori A and 
1 Ori C hardly comparable. The detected magnetic field 
is indeed much weaker than the critical limit above which 
MHD instabilities are inhibited (about six tim es the equipar- 
tition field or 600 G in the case of £ Ori A, lAuriere et al.l 
l2007f > and may thus result from exotic dynamo action; the 



non-dipolar nature of the detected field could be additional 
evidence in favour of this interpretation, fossil fields being 
expected to have very simple topologies in evolved stars. Ad- 
ditional spectropolarimetric observations of £ Ori A at dif- 
ferent epochs (searching for potential variability of the large- 
scale field) and of similar 'normal' rapidly-rotating stars are 
obviously necessary to explore this issue in more details. 

C omputing the wind magnet ic confinement parameter 
7/* of lud-Doula fc OwockH |2002l 'l for £ Ori A and taking 
B ~ 30 — 50 G (at the magnetic equator), R — 25 Rq, 
M = 2 x 10~ 6 M yr _1 and v x = 2,100 kms" 1 (see 
Sees. [3] and g} yields 77, ~ 0.03 - 0.07. The magnetic field of 
£ Ori A is therefore just strong enough (according to theo- 
retical predictions) to start distorting the wind significantly 
|ud-Doula fc Owockil liobl ). The observed rotational mod- 
ulation in Ha, H/3 and the C ill lines confirms this first 
conclusion; the variation in mass loss rate that we measure, 
corresponding to a density contrast of ~1.4, is compatible 
with what numerical simulatio ns of magnetically confine d 
winds predict (see, e.g., Fig 8 of lud-Doula fc Ow ocki 2002). 

Note also that the observed line blueshift and asymme- 
tries (cf. Sec. I3.2l and Fig. [7] for an illustration on the case of 
He 1 492 nm) are (rotation) phase dependant. The observed 
blueshift of most lines is maximum when the magnetic poles 
(i.e. the open field lines) cross the line of sight (at phase 
0.8 and 0.45, see Fig. |7J; more data (collected in particu- 
lar over a longer baseline and densely sampling the rotation 
cycle) are of course needed to confirm this and to specify 
how exactly the line shifts and shape are varying with rota- 
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tion phase (e.g. with 2 unevenly spaced maxima in the line 
blueshift per rotation period). It however suggests already 
(?) that the line profile variations reflect the varying condi- 
tions in which the wind form at the surface of £ Ori A (as a 
result of the varying local field topology over the star) and 
(ii) that these variations can potentially be used to trace 
the density at the base of the wind and its variations with 
the local magnetic topology over the surface of the star. 

On both 9 1 Ori C and HD 191612, Ha is rotationally 
modulated as a result of the magnetic obliquity (with respect 
to the rotation axis) , with maximum emission occuring when 
the magnetic pole comes closest to the observer. Similarly, 
maximum absorption in UV lines (with highest blueshifted 
velocities) are observed when the magnetic equator is cross- 
ing the line of sight. Extrapolating these results to £ Ori A, 
we would have first expected Ha in £ Ori A to show maxi- 
mum emission twice per rotation cycle, at phases of about 
0.40 and 0.90 (see Fig. |9}, in contradiction with what we 
observe; while Balmer emission indeed peaks at phase 0.95, 
phase 0.51 rather corresponds to minimum (rather than 
maximum) emission (see Fig. [7}. 

The analogy with 8 1 Ori C and HD 191612 can obvi- 
ously not be directly applied to £ Ori A. Given the much 
weaker wind magnetic confinement parameter of £ Ori A 
(roughly equal to 10 for both 6 1 Ori C and HD 191612), 
this is not altogether very surprising. In particular, the 
Alfven radius is much closer to the surface of the star in 
£ Ori A, probably not further than 0.05-0.1 -R* above the 
surfacqj (as oppos ed to 1 R+ above the surface for 8 1 Ori C, 
iDonati et al.ll2002f). In the magnet i cally confined wind- shock 
model (|Babel fe Montmerld 1 19971 ; IDonati et aT] l2002h . the 
rotational modulation of Ha emission can be ascribed, in a 
generic way, to the varying aspect of the magnetic equato- 
rial disc up to the Alfven radius; in the case of £ Ori A, this 
variation is expected to be minimal. We speculate that most 
of the redshifted Ha emission comes from a region located 
just above the photosphere (at the very base of the wind) 
and essentially reflects a difference between both magnetic 
poles (strongest emission being observed in conjunction with 
the weakest magnetic pole, see Fig. [5] bottom panel). 

We also note that the excess absorption that both 
Balmer lines exhibit twice per rotation period in their dis- 
tant blue wing (see Fig. [7] left and right panels) behave as 
UV absorption lines do in 9 1 Ori C; we propose that they 
reflect the magnetic equator crossing the line of sight (at ro- 
tation phases of about 0.20 and 0.75). The maximum radial 
velocities associated to these absorption components (up to 
about 500 kms -1 , i.e., less than 0.25«oo) confirm that they 
correspond to material located within the Alfven radius. 
More data (densely sampled over several rotation cycles) 
are needed to investigate this issue more closely, and to pin 
down unambiguously the origin of the various Ha and H/3 
components. 

The 569.6 nm C ill double-peak emission line is 
also a significant difference with respect to 8 1 Ori C and 
HD 191612 (where the line only features a single peak emis- 



2 The corotation radius, i.e., the radius at which the Keplerian 
period equals the rotation period at the surface of the star, is 
equal to about 2 i?* in ( Ori A, i.e., 1 R* above the surface of 
the star. 



sion). The observed modulation is apparently related to 
the magnetic topology, with the red emission peaking at 
phases of maximum magnetic field and the blue emission 
peaking when the magnetic equator is crossing the line of 
sight (both phenomena occuring twice per rotation period). 
The maximum velocities of both components (up to about 
200 kms" 1 ) also argue for the formation of this line within 
the Alfven radius and, therefore, it is likely a good indica- 
tor of the influence of the magnetic field on the launching 
of the wind. Further observational and theoretical studies 
are again required to examine how this line responds to a 
magnetised wind. 

At the very least, our results demonstrate that the mag- 
netic field of £ Ori A has a significant impact on the wind 
despite being below pressure-equipartition and the weakest 
detected ever in a hot star. Given the obvious importance of 
this result for our understanding of massive magnetic stars, 
we need to confirm and expand the present analysis with 
new data collected over several rotation periods of £ Ori A, 
i.e., over a minimum of 20 nights; renewed observations will 
indeed allow us (i) to obtain an accurate measurement of the 
rotation period, (ii) to derive a fully reliable modeling of the 
large-scale magnetic topology and (iii) to estimate whether 
the field is intrinsically variable as usual for dynamo topolo- 
gies, e.g., on a typical timescale of 1 yr, and (iv) a detailed 
account of how wind lines (and in particular Ha, H/3 and 
the 569.6 nm C ill lines) are modulated with the viewing 
aspect of the magnetic topology. 
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